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Abstract: We complete the ‘paraxial’ (small-angle) ray optics cloaking 

formalism presented previously [Choi and Howell, Opt. Express 22, 29465 
(2014)], by extending it to the full-field of light. Omnidirectionality is then 
the only relaxed parameter of what may be considered an ideal, broadband, 
field cloak. We show that an isotropic plate of uniform thickness, with 
appropriately designed refractive index and dispersion, can match the phase 
over the whole visible spectrum. Our results support the fundamental limits 
on cloaking for broadband vs. omnidirectionality, and provide insights into 
when anisotropy may be required. 
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1. Introduction 

Invisibility efforts by scientists have focused on building a broadband, omnidirectional cloak¬ 
ing device HI. Such a cloak that works in macroscopic dimensions, for the entire visible spec¬ 
trum, with the full-field of light may be considered an ‘ideal’ cloaking device ||2l. Much of 
this research was spawned by the initial omnidirectional cloaking designs, which used artifi¬ 
cial materials (called ‘metamaterials’) to control electromagnetic fields/waves Ell. This field 
was called ‘transformation optics,’ since coordinate transformations were used to design the 
necessary material properties that morphed space for light. Transformation optics has produced 
many ingenious designs and devices 0. However, creating the required spatial distribution of 
metamaterials, their anisotropy, and causality, limit the spectrum to single frequencies, or to a 
narrow bandwidth, and is particularly difficult to manufacture for visible frequencies 010. 

Much of the scientific work on cloaking has retained the omnidirectionality and full-field 
(amplitude and phase of light waves) nature of transformation optics 121 . With this, to make 
forays into large bandwidths in the visible spectrum, researchers have built reflecting ‘carpet 
cloaks’ Em, proposed non-Euclidean geometry mapping 0, demonstrated macroscopic ray 
optics cloaks (matching directions, rather than phases, of fields) Em, and suggested strongly 
diamagnetic superconductors HI. 



Simultaneously, researchers have investigated the fundamental limits on bandwidths for these 
omnidirectional cloaks. It is generally agreed that causality requires such cloaking of a non-zero 
volume, either passive or active ifT^ . to be only possible for single frequencies iflBll . Chen et 
al. showed that zero scattering cross-section of a cloaked object (for an omnidirectional trans¬ 
formation optics cloak in two-dimensions (2D)) cannot be attained by neighboring frequencies 
of the target frequency, no matter how narrow the bandwidth m. They then derived an upper 
bound for the bandwidth that was proportional to the scattering cross-section radius. Hashemi 
et al. have also provided a theoretical bound for transformation optics cloaks, even if perfectly 
manufactured, which shows increased imperfections for increasing bandwidth ca. Derived 
from causality, their “diameter-bandwidth product” limit states that the effective bandwidth 
is inversely proportional to the cloaked object diameter. They suspect a similar sensitivity to 
‘carpet cloaks’ as well. Similarly, Monticone and Alii showed that any linear, passive, non- 
diamagnetic cloak must increase the overall scattering, compared to an uncloaked object, when 
integrated over all wavelengths ifThl . Even if scattering is suppressed in a finite bandwidth, this 
will need to be “paid back” in the rest of the spectrum. On the other hand, Leonhardt and Tyc’s 
broadband solution is a three-dimensional (3D), anisotropic cloak 0. Their “price to pay for 
practical invisibility” is relaxing the full-held requirement to allow time delays. 

Our approach has been to begin with broadband cloaking for the entire visible spectrum, 
but within the paraxial approximation, then to expand the cloaked size and angles ifTTl . This 
is typically the method used in optical engineering, as increasing the held-of-view (angles) 
and/or numerical aperture (refractive index multiplied by the solid angle of the accepted light) 
can be challenging M- So in this paper, we use the dehnition for “cloak” that is to “hide,” 
since the cloaked object is not necessarily covered in all directions as if by a garment. Our 
previous ray optics demonstration relaxed both angles and phase to show that only four lenses 
were needed for 3D, macroscopically scalable, stand-alone cloaking for small angles OH. Our 
paraxial cloak was the simplest design we found that used isotropic optics designed for the 
entire visible spectrum. These optics have been well understood for centuries ll20ll and are 
readily available off-the-shelf. Interestingly, Hashemi et al. concluded by suggesting a relaxing 
of angles and phase, to possibly obtain a precise understanding of when the cloaking problem 
becomes easy ca. They considered knowing what cloak was the “weakest” relaxed, yet still 
practical for large objects, to be valuable, due to the difficulties of transformation optics cloaks. 

2. Extending paraxial cloaking to include the full-field 

We now remove the ray optics requirement, to complete a paraxial cloaking theory that works 
for the full-held (matched amplitude and phase). Such paraxial full-held cloaking satishes all 
but the omnidirectionality condition of an ‘ideal’ cloaking device. This ‘paraxial’ formalism 
can work for up to ±30° ED, which is practical for many cases where the cloak is not placed 
immediately before an observer. Although we discuss the propagation of a monochromatic 
held, since an arbitrary held of light can be written as a linear superposition of monochromatic 
waves, our theory extends to broadband without loss of generality. 

Duan et al. recently provided phase matching with their unidirectional cloaking system based 
off of geometric optics ll22l . They provided a heuristic reason why rays that pass through the 
edge of their optical system may match the phase of rays that pass through the center instead, 
for their particular setup with split lenses. They then simulated phase-matched cloaking for 
discretely separated wavelengths of micro-waves to mm-waves. Here we analytically show 
phase matching for general paraxial optical systems, including continuously multidirectional 
cloaks. It is based on the formula given by Siegman and others, for propagation of any paraxial 
field through a generalized paraxial optical system ED- The resulting formula is proved using 
Fermat’s principle and Huygens’ integral, by accounting for the optical path lengths of all rays. 



Fig. 1. Ideal (‘perfect’) paraxial cloak. The image from a cloaking device is the same as 
the object. Propagation of light rays can be described by an ‘ABCD’ matrix, n is the index 
of refraction of the ambient medium. L is the longitudinal length of the device along the 
center z-axis. £’i,i?2 ars ths input and output fields, respectively, at z = zi and Z 2 - Phase 
matching is achieved with a flat plate with index and length L^. 


Field propagation based off of Huygens’ principle of wavelet propagation, is effective and 
widely used in diffraction theory and Fourier optics ||2^ . We first assume that the ambient 
medium is spatially uniform with index of refraction n, and that the optical system is rota- 
tionally symmetric for simplicity. Non-uniform medium (US) or non-rotationally symmetric 
systems m can be derived from here. 

Additionally assuming no limiting apertures, Huygens’ integral in the Fresnel (or, paraxial) 
approximation is given by lElEl (SeeFig.IB: 

~ f ~ f TT 0 0 0 0 1 

E 2 {x 2 ,y 2 )= J j Ei{xuyi)exp^i—[A{xi+yi)-2{xiX2+yiy2)+D{xj+yj)]\dxidyi. 

( 1 ) 

A,B,C,D are the ABCD matrix coefficients. Lq = Y.i’^iEi is the on-axis optical path length, 
where each optical element has index of refraction «,■ and physical thickness L,- along the 
longitudinal axis (z). Aq and Aq are the free space wavelength and wave vector, respectively. 
El ,£2 are the complex, spatial amplitudes of the input and output field distributions, respec¬ 
tively {Ei,E 2 in Fig.[T]are their real parts, but without the harmonic time dependence). 

In our previous work, we stated that a ‘perfect’ cloaking device (of length L) simply replicates 
the ambient medium throughout its volume, and that its ABCD matrix is given by mi 
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Thus, the propagated field for a ‘perfect’ full-field cloak, of length L, is: 

E^°^{x2,y2)= '^^^^ j j Ei[xuyi)e^^iyi^[{x\+y\)-2{xiX2+y\y2)+{xl+yl)^Axidyi. 

(3) 

A ray optics cloak satisfies Eq. (I2]i. By comparing Eq. ([T]) and Eq. Q, we see that a ray 
optics cloak can be a full-field cloak, if e^'^ 0^0 = g-<konL^ Specifically, this is the case when 
koLo = konL (absolute phase-matching), or when Lq = nL (mod Ao) (phase-matching to integer 
multiples of 2;r). The significance of these conditions is that they allow for phase-matched, full- 
field cloaking for any incoming fields within the paraxial approximation. The phase-matching 
condition only needs to be satisfied once for a given optical system, and then all ray directions 
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Fig. 2. Dispersion of phase-correcting plate for a four lens symmetric cloak. Various 
thin, flat plates are used to match the full-field cloak condition (Eq. (O), based on the ray 
optics four lens cloak in OH- The dispersion relation of the refractive indices for these 
plates are shown, for various values of m and L^, using Eq. 0 in ambient air. 

and positions, or any field distribution (Ej) will exit as if traversed through ambient space. Note 
that a ray optics cloak will usually be phase-matched with a full-held cloak, by integer multiples 
of 27r, for multiple, but discretely separated, wavelengths automatically. However, we can do 
better and match it for a continuous, broad bandwidth with appropriate dispersion control. 

One method for broadband phase matching is to add a thin, flat plate to a ray optics cloak, 
anywhere between the background object(s) and the observer. Since it is “thin” and flat, the 
original ABCD matrix (Eq. (|2]i) will be unchanged ll25l l26l and only the factor outside 

the integral of Eq. O will be affected. We consider the case shown in Eig. [T] where the flat 
plate is placed immediately after a ray optics cloak. Let L' = be the total length of the 

original ray optics cloak (so L = Lc+L'), where N is the number of original optical elements. 
So Lo ncLc + Y!i=\ in Eq. ([T]), and nL n{Lc + L') in Eq. Q. Here, are the index of 
refraction and physical longitudinal length, respectively, for the optical element, and iicLc 
are likewise for the correcting thin plate. It can be shown that a plate is “thin” when \Lc/L'\ <C 1 
and \Lc/L'\ ^ |«c/k|- For our thin correcting plate, its phase-matched dispersion relation is 



(4) 


m is the integer multiple of 2;r for which the phase is matched, with discrete but infinite choices. 

There are infinitely many combinations of Lc and m that are possible for nc{Xo,m,Lc). Ul¬ 
timately, the choice of which to use should depend on physically realizable ric’s. A few such 
solutions, that phase-corrects the four lens ray optics cloak in lfT9l . are shown in Eig. |2] Eor 
simplicity, we assumed that the cloaking system is placed in air, with n = = 1 for all wave¬ 

lengths, and we ignored any coatings. The lens coatings can be modeled as additional thin 
plates as needed. We also only consider the visible spectrum (400-700 nm, for our discussion). 
A thicker compensating flat plate can reduce the required dispersion range, but this also affects 
the imaging quality, so we have purposely limited its thickness here. 

Many of the solutions for positive refractive indices require anomalous dispersion, as shown 
in Eig. 12] Eor the purposes of cloaking, such dispersion for broadband spectrum with low losses 
is needed to imitate ambient space properly. This is similar to the finite bandwidth cloak made 
of anisotropic layers, as suggested by Kildishev et al., which required strong anomalous dis¬ 
persion combined with loss compensation EH. Costa and Silveirinha have suggested using 
nanowire metamaterial to provide such anomalous dispersion, and they numerically calculated 







their index of refraction to be near 2.9 ~ 3.1 for the entire visible spectrum ||28l. They can 
achieve these low loss, broadband, anomalous dispersion properties by utilizing the collective, 
spatial properties of metamaterials. This allows the high loss and narrow band properties of 
typical transparent materials, imposed by the Kramers-Kronig relations, to be overcome 12^ . 
In addition, Theisen and Brown have experimentally demonstrated anomalous dispersion for 
0.5 ^ 1 jim wavelengths, with Gallium implanted Silicon pads 1^ . These have refractive in¬ 
dices near 2 or 4, with variations of about 1 over the spectrum, depending on the doping level. 

Negative-index metamaterials are good candidates for phase-correction as well ED. The uti¬ 
lization of metamaterials may have been expected, since we expanded ray optics cloaking to 
the field cloaking aspect of transformation optics. Much progress is being made, both theoret¬ 
ically and experimentally, for creating negative index materials for broadband optical frequen¬ 
cies EH. Some refractive index values demonstrated experimentally include between -3 and 0 
for microwave frequencies in 2D ES, n = 0.63 at 1,200 nm to n = -1.23 at 1,775 nm in a low 
loss, 3D bulk material ED . and n = 1 to -7.5 for 1.1-2.4 fim wavelengths ESl . 

3. Discussion 

By relaxing only omnidirectionality for an ‘ideal’ cloak, we have shown how to match the phase 
for the whole visible spectrum. The phase-matching plate may require negative index metama¬ 
terials or anomalous dispersion, which are broadband and low loss, but current research has 
shown much progress in this regard. We had shown that building a 3D, broadband, macroscopic 
cloak, for the visible spectrum, can be fairly “easy” for ray optics in the small angle limit US). 
Extending this to the full-field seems to not be too challenging with recently developed mate¬ 
rials. As with typical lens designs, we expect that extending such cloaks to large angles may 
be difficult, though possible. However, making it work for all angles appears to likely be fun¬ 
damentally limited. Realizing omnidirectionality from a paraxial full-held cloak has not been 
achieved to our knowledge, as this coincides with achieving broadband with transformation 
optics cloaks. By showing paraxial cloaking that is broadband to be practical, but without a 
similar ease for omnidirectional and broadband cloaks, our work supports recent work showing 
a trade-off between broadband and scattering cross-section. This is because small scattering 
cross-section implies large angles, and zero cross-section gives omnidirectionality. 

It is interesting that paraxial full-held cloaking does not require anisotropy, though being 
3D, macroscopic, and broadband. Anisotropy seems to be a requirement when creating an om¬ 
nidirectional cloak, from a paraxial cloak, not necessarily a property of 3D or held cloaks 
alone El HI El [361. Although our cloak is broadband, Greenleaf et al. have shown that the 
cloaked material itself has eigenfrequencies that need to be avoided for proper cloaking ETlESlI . 
As a side, most, if not all, cloaks to date, can be detected by measuring a pulse incident on the 
cloak EllH ■ This is the same for our cloak presented here, unless absolute phase is matched with 
negative index materials. Finally, an open question seems to be whether an isotropic, 3D, om¬ 
nidirectional, broadband cloak can be achieved for ray optics. This is possible for anisotropic 
materials El. but an isotropic, 3D, omnidirectional, full-held cloak is not possible E9l l40l . 
However, accepting time delays might allow some freedom, as shown here where ray optics 
allowed relaxing of the material requirements. 
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